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The phenomenon of shape coexistence is discussed within the self-consistent Hartree-Fock method 
and the nuclear shell model. The occurrence of the coexisting configurations with different intrinsic 
shapes is traced back to the properties of the effective Hamiltonian. 
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I. INTRODUCTION 

The phenomenon of nuclear coexistence manifests it- 
self in the presence of close-lying nuclear states with very 
different intrinsic properties. Spectacular examples of co- 
existence are superdeformed states, low-lying deformed 
states in spherical nuclei, high-_R' isomers, and pairing 
isomers. One of the most exciting aspects of the coexis- 
tence phenomenon is the fact that the coexisting excited 
states often retain their identity at rather high excitation 
energies. Indeed, the lowest observed gamma transitions 
in superdeformed bands correspond to states lying sev- 
eral MeV above the yrast line, i.e., in the region of very 
high level density. This means that superdeformed con- 
figurations do not mix with many near-lying states; they 
are very diabatic. 

In most cases, energies of coexisting states strongly 
depend on particle numbers. For instance, the neutron- 
deficient Hg isotopes have well-deformed prolate ground 
states containing the high-j proton orbitals "intruding" 
across the Z—82 gap while the ground states of heavier 
Hg isotopes are only weakly deformed and they can be 
associated with oblate shapes. Such crossings between 
coexisting structures are particularly interesting in light 
nuclei; they can give rise to the presence of deformed 
ground states in magic nuclei such as ^^Mg. For an ex- 
tensive review of shape coexistence in light and heavy 
nuclei, we would like to refer the reader to Refs. 

In the mean-field approach, the mean fields, in which 
nucleons move as independent (quasi)particles, can be 
obtained from a knowledge of the effective forces acting 
between nucleons using the Hartree-Fock (HF) theory. 
For a particular choice of nucleon-nucleon force and pro- 
ton and neutron numbers, the mean field may be spheri- 
cal or deformed. The minimization of the HF energy un- 
der the constraint of holding any given nuclear multipole 
moment fixed can be carried out over a range of collective 
parameters. When more than one local minimum occurs 
for the total energy as a function of deformation, shape 
coexistence may result. 

The phenomenon of shape coexistence can be traced 



back to the nuclear shell structure. First, the sizes of 
spherical magic and semi-magic gaps in the energy spec- 
trum determine the relative positions of many-particle 
many-hole intruder configurations with respect to the 
ground state. Secondly, the strength of the particle- 
vibration coupling responsible for the breaking of spher- 
ical symmetry (i.e., the development of deformation) is 
strongly dependent on the relative distance between indi- 
vidual shells 1^,^ . This can be qualitatively understood 
by means of Strutinski energy theorem j5|-|^ , which states 
that the total HF energy can be written as 



£^shcn + 0('5p'), 



where E is the smoothed HF energy. 
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is the shell energy (shell correction), h^^ is the HF 
Hamiltonian, p is the single-particle density, and Sp rep- 
resents the contribution to p due to shell effects. The shell 
energy reflects single-particle properties of the Hamilto- 
nian; it lowers the binding energy if the Fermi level is 
situated in the region of low level density |^. Hence, 
-BshcU is very sensitive to the details of the single-particle 
spectrum. The average HF energy E behaves, roughly, 
like that of the nuclear liquid drop ||] . In contrast to the 
shell correction term, it reflects the average properties of 
the interaction. In the geometric picture, deformation 
comes from the subtle interplay between E and iSshoii- 

It is important to realize that, in the mean field ap- 
proach, the excitation energy of the coexisting (intruder) 
state involves a difference between the binding energies 
of coexisting minima, hence it can easily be obscured 
by a different quality of the theoretical description for 
these states. For instance, in recent work |^], the self- 
consistent mean-field theory was used to explain at the 
same time the two-particle separation energies in the first 
and second wells, and the excitation energies of superde- 
formed states. While for the particle separation energies 
in the first and second wells a good agreement with ex- 
perimental data could be found, this did not hold for 
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relative energy differences between the wells. This exam- 
ple demonstrates that while the intrinsic configurations of 
coexisting states do not depend in a very sensitive way on 
details of calculations, the relative energies of coexisting 
states are strongly affected by model uncertainties such 
as treatment of pairing, surface tension, or level density, 
i.e., those particular properties of force parameterizations 
that determine the deformability of a nucleus. 

Another source of uncertainty in HF calculations con- 
cerns the corrections which should be added to the calcu- 
lated energies to account for dynamical correlations asso- 
ciated with zero-point fluctuations. Since the dynamical 
corrections can be different in coexisting minima, they 
can influence the predicted excitation energies Q . 

Thus far, we have discussed shape coexistence from 
the viewpoint of the mean-fleld theory. An alternative, 
and in many respects complementary, approach is pro- 
vided by the nuclear shell-model calculations which aim 
to determine the fully correlated states in an expansion 
basis of the few active shells near the Fermi energy. Here 
also, the interplay of the shell structure in nuclei and 
the deformation-driving proton-neutron residual interac- 
tion is the key to understanding the shape coexistence 
in terms of the spherical shell model. The main mecha- 
nism for shape coexistence here is the multiple particle- 
hole (or pair) excitation across the closed shell By 
promoting nucleons out of the closed core into the next 
higher shell, nucleons of both kinds interact to produce 
deformed structures. This approach to shape coexistence 
was first invoked to explain the 0"*" state at 6.05 MeV 
in 1^0. 

In the shell-model framework, the band-head of the 
deformed intruder configuration can be written as : 



E* = AE, 
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where the contributions to the intruding coexisting state 
energy are: AEp\i ~ the unperturbed particle-hole exci- 
tation energy; Ai?pair - the change in pairing correla- 
tion energy resulting from the particle-hole excitation; 
AEm - the change in (proton-neutron) monopole in- 
teraction energy; and AEq ~ the change in (proton- 
neutron) quadrupole interaction energy. In Ref. [Q, 
the AE-pYi and ASpair terms were estimated using ex- 
perimental separation energies, and the AEm and AEq 
terms were calculated using a delta interaction for AEm 
and the schematic quadrupol e-qu adrupole interaction for 
AEq. When applying Eq. (1.3) to pair excitations in 
closed-shell nuclei, it was found that the pairing en- 
ergy, monopole energy, and quadrupole interaction en- 
ergy tend to reduce Ai?ph. For a fixed value of N, the 
monopole energy monotonically decreases with Z, and 
the quadrupole term reaches its maximum at midshell. 

The main objective of this paper is to trace the phe- 
nomenon of shape coexistence to properties of the ef- 
fective nuclear Hamiltonian. In our analysis, we apply 
the mean-field self-consistent Skyrme-HF method and 
the nuclear shell model. The sensitivity of the inter- 
play between the coexisting configurations is discussed in 



terms of several key quantities such as the single-particle 
splitting, pairing correlations, and surface tension (in 
the mean-fleld approach), as well as the monopole en- 
ergy, quadrupole correlation energy, and effective single- 
particle spectra (in the shell model). 

The paper is organized as follows. The four regions of 
shape coexistence studied in this work, namely the de- 
formed iV=20 and A^=28 regions, the region around *°Zr, 
and the region around ^^Zr, are briefly reviewed in Sec. ||. 
The results of the mean-field and shell-model analysis are 
contained in Sees. |l| and respectively. Finally, the 
main conclusions of this work are summarized in Sec. 



II. EXAMPLES OF SHAPE COEXISTENCE 
A. The deformed A'^~20 region 

The neutron-rich nuclei with A^«20 are spectacular ex- 
amples of coexistence between spherical and deformed 
configurations in the sd shell (8< Z,N < 20). A classic 
example is the magic nucleus nnMg, which has a very low- 
lying 2+ state at 886 keV [|5[ and an anomalously high 
value of the two- neutron separation energy S2n- Ground- 
state deformations in this mass region can also explain an 
anomalous isotope shift in '^^N a [17| and a major decrease 
in S2n in 3i.33,35]v^g^ g^j^^ 30]v^g fl^^ . The large deforma- 
tion of ^^Mg has been inferred from the intermediate- 
energy Coulomb excitation studies [Q. 

In many calculations based on the mean-field theory, 
deformed ground states have been predicted in nuclei 
from the ^^Mg region (sometimes dubbed as an "island 
of inversion"). In the early Skyrme-HF calculations of 
Ref. ijl^ with the SIII and SIV interactions, large prolate 
deformations in ^^'^^Na were obtained and explained in 
terms of neutron excitations from the sd shell to the /y /2 
shell. Similarly, low-lying 2+ states have been predicted 
in A^=20 nuclei based on the energy density formalism 
p9[ . A sudden onset of large ground-state deformations 
(^~0.3-0.4) around iV=20 was also predicted in the cal- 
culations based on the macroscopic-microscopic method 

013 ■ 

In the shell-model language, the structural changes 
around '^^Mg can be attributed to the cross-shell particle- 
hole excitations to the /7/2 shell. Early shell-model calcu- 
lations in a rather restricted configuration space (no more 
than two neutrons in the /7/2 shell) |22| were able to re- 
produce the increased quadrupole collectivity at N=20. 
A similar conclusion was drawn in other shell-model cal- 
culations in the {sdpf) model space [p3|-p7[, allowing only 
two-neutron particle-hole excitations from the sd to the 
pf shell , an d in the schematic analysis of Ref. based 
on Eq. (Ilj). 

Recently, the onset of deformations in this region has 
been a subject of much theoretical work, strongly moti- 
vated by the prospects of detailed experimental spectro- 
scopic studies at ISOLDE In Ref. |30| (see also Ref. 
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based on the RMF theory with the NL-SH force and 
the constant gap BCS treatment of pairing, "^^Mg was 
calculated to be spherical. They concluded that this re- 
sult, together with the previous RMF study with the NLl 
force P4], did not depend on the choice of the RMF pa- 
rameterization. Spherical ground-state deformation for 
■^^Mg has also been obtained in the Skyrme-HFB cal- 
culations with the SIII, SLy4, and SkP forces and with 
a density-dependent zero-range pairing interaction [ p3[ , 
in the HFB study based on a Brueckner G-matrix de- 
rived from a meson-exchange potential with the density- 
dependent meson masses [p4[ , and also in the Gogny- 
HFB calculations of Ref. |3| . The authors of Ref. Q 
noted, however, that strong deformation effects around 
■^^Mg could appear due to dynamical correlations. Their 
collective wave functions of -^^Ne and ■^■^Mg, calculated 
with the collective Hamiltonian, have pronounced max- 
ima at large deformations. 

B. The deformed iVf=s28 region 



C. The Af«Z«40 region 

The proton-rich N=Z nucleus ^'^Zi lies in the center 
of the well-deformed A«80 region |^^. The sizable en- 
ergy gap at particle number 40 separates the pf spher- 
ical shell from the (79/2 orbital. However, this spherical 
subshell closure is not sufficiently large to stabilize the 
spherical shape. Experimentally ^"Zr seems to be a 
well-deformed rotor. According to the mean-field theory, 
this is due to the presence of the deformed single-particle 
gap at N,Z=40- the resulting deformed shell effect turns 
out to be stronger than that at the spherical shape. 

Microscopic calculations based on the symmetry- 
projected variational model [Q, Skyrme-HF theory 
j|9|,|o|, RMF theory and the restricted-space HFB 
calculations predict a deformed ground-state mini- 
mum for ^^Zr. Only in a very few calculations, such as 
the RMF calculations with the NLl parameter set [ p3[ , 
was a spherical ground state obtained. (See, however, 
the discussion in Ref. jsill .) 



Another, recently discovered, island of inversion are 
the neutron-rich nuclei from the pf shell centered around 
^'*S28. Experimentally, /3-decay properties of ^"^S and 
45-47q]^ have been studied in Refs. [^,0. Based on the 
QRPA analysis of measured half-lives, it was concluded 
that **'*S was deformed. This has been confirmed recently 
in a series of intermediate-energy Coulomb excitation 
studies |]38| , |39| which revealed rather large B{E2; 0+^ 
2i) values in this neutron-rich region, suggesting a sig- 
nificant breaking of the A^=28 core. (For recent mass 
measurements around "''^S, see Ref. pO[|.) 

The HF+BCS calculations with the Skyrme interac- 
tions SHI and SkM* and the RMF calculations with 
the parameter set NL-SH |^,^ predicted the apprecia- 
ble breaking of the N=28 core and deformation effects 
around '''*S. A later study Q, based on the RMF ap- 
proach without pairing and using the TMl parameter set, 
predicted the neutron-rich sulfur isotopes to be deformed. 
In particular, '*'*S was found to be prolate, in agreement 
with the RMF results of Refs. QQ. Only very re- 
cently, the interplay between deformed mean-field and 
pairing correlations in this mass region has been prop- 
erly considered in the framework of relativistic Hartree- 
Bogolyubov (RHB) theory [4^ using the NL3 effective 
interaction for the mean-field Lagrangian and the Gogny 
interaction DIS in the pairing channel. Again, deformed 
shapes around "''^S have been calculated. 

An erosion of the N—28 gap in the sulfur isotopes 
has also been found in shell- model calculations |45| per- 
formed in a large configuration space (the full sd shell 
for protons and the full pf shell for neutrons). The au- 
thors concluded, however, that the shell-breaking effects 
around ^"^S were much weaker compared to the NRi20 
neutron-rich region. 



D. The N^56, Z^40 region 

Nuclei from the heavy-Zr region (Z?^40, N>56) exhibit 
a wealth of coexistence phenomena The strong 

dependence of observed spectroscopic properties on the 
number of protons and neutrons makes the neutron-rich 
AwlOO nuclei a very good region for testing various mod- 
els. Theoretically, strong shape variations in this region 
may be attributed to shell effects associated with large 
spherical and deformed subshell closures in the single- 
particle spectrum js^ ]. 

According to calculations based on the mean-field ap- 
proach, the occupation of the /111/2 neutron and 179/2 
proton orbitals is essential for understanding the de- 
formed configurations near ^^'^Zr [s^ . The best exam- 
ples of shape coexistence in this region are the Sr, Zr, 
and Mo isotopes with iV«58. In the language of the 
deformed shell model, the onset of deformation around 
N—58 can be associated with the competition between 
the spherical gaps at Z=38, 40, and A^=56, and the de- 
formed subshell closures at particle numbers Z=38, 40, 
and iV=60, 62, and 64. Theoretically, the deficate energy 
balance between spherical and deformed configurations 
depends crucially on the size of these gaps. As discussed 
in Refs. [pl Js^JS^ , the deformation onset at NKi58 re- 
sults from the subtle interplay between the deformation- 
driving neutron-proton quadrupole interaction and the 
symmetry-restoring monopole force responsible for shell 
effects. 

Equilibrium deformations and moments, potential en- 
ergy surfaces, the microscopic structure of coexisting con- 
figurations, and shape transitions in the heavy-Zr region 
have been calculated by maiw authors. (For an extensive 
list of references, see Ref. [p5|.) In most cases, the cal- 
culations show large deformations in the Sr, Zr, and Mo 
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isotopes with A''>60. The details of the shape transition 
near N=58 are, however, predicted differently by various 
models, the onset and rapidity of this transition being 
very sensitive to the actual parameterization used |49|,p| . 
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-61 + 62 - K + 6^2 J 



6i + 2fo2-^fo'i + 36^Wj 



(3.6) 



III. SKYRME-HARTREE-FOCK CALCULATIONS 

A. The Skyrme Energy Functional 

Our implementation of Skyrme forces is based on the 
standard ansatz as it has now been used for more than 
two decades The total binding energy of a nucleus 
is obtained self-consistently from the energy functional: 
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£c = I d?rd\'p,{r)^-^p,{r'), (3.4) 



£sk,is is the spin-orbit functional, fpair is the pairing en- 
ergy, and f CM is the center-of-mass correction. 

The functional employs the usual particle densi- 
ties pq = X^QGg "'"l^"!^' kinetic densities Tq = 
X^qgij "-ctl^^ctP' ^^"^ the spin-orbit densities Jq = 
X^QGg^^^a"" ^ ^V'a, where ipa are the single-particle 
(canonical) wave functions and q stands for either pro- 
tons or neutrons. The total isoscalar density is p=Pp+Pn 
and similarly for r and J . The ria—v'^ is the BCS occu- 
pation weight (see below). 

The terms discussed above are always defined in the 
same way for all Skyrme parameterizations. This is 
not the case for the remaining terms in Eq. ( |3.l|) . The 
Coulomb-exchange functional is usually treated in the 
Slater approximation 
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(3.5) 



but it is omitted in definitions of some published Skyrme 
forces. All the parameterizations considered in this work 
require this term. 

The spin-orbit functional can be written as 



■'Sk,ls 



^64/0 V • J 
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This spin-orbit functional encompasses two different op- 
tions, namely, one either ignores the contributions 
(6*18 = 0) or takes them into accou nt (0\s = 1). Fur- 
thermore, the spin-orbit functional (3.6) is given in the 
extended form of |^ which allows a separate adjust- 
ment of isoscalar and isovector spin-orbit force. The stan- 
dard Skyrme forces use the particular combination 64=64 
which was motivated by the derivation from a two-body 
zero-range spin-orbit interaction |Q, but these partic- 
ular settings are not obligatory when taking the view- 
point of an energy-density functional. Thus, various op- 
tions exist in the published literature, and we shall use all 
combinations of them in the examples discussed in this 
work. 

Similarly, there are basically two different options for 
handling the center-of-mass correction, i.e.. 



' 2m 



ft 
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(3.7) 



For dcM—Q, the center-of-mass correction is implemented 
before variation by the simple trick to let the nucleon 
mass m — > mjod = m — m/A. The option 0cm=1 uses a 
more correct expression, but it is difficult to implement 
in the fully self-consistent manner due to the two-body 
nature of the P^m operator. Hence standard parame- 
terizations with 9cM—^ apply this correction after vari- 
ation for the given mean-field solutions obtained with 
the center-of-mass correction ignored. (For examples of 
a fully self-consistent treatment, see Ref. [^.) 

Since there are more than 80 different Skyrme param- 
eterizations on the market, the question arises, which 
forces should actually be used when making predictions 
and comparing with the data? An extensive list of forces, 
together with their properties, can be found in Ref. [|63| . 
To select a manageable number, we have computed the 
overall quality factor which reflects the predictive power 
of the force for the basic ground-state properties (masses, 
radii, surface thicknesses), and confined further analysis 
to the best performing parameterizations. From these, 
we have chosen a still smaller subset with sufficiently dif- 
ferent properties to explore the possible variations among 
parameterizations. This subset contains: SkM* [ |64[ , 
SkT6 |6|, lei, SkP SLy4 and Skll, SkI3, 
and SkM from Ref. We have also added two addi- 

tional forces from a recent exploration [Q. These two 
forces are labeled SkO and SkO'. A list of the parameters 
for these forces is given in Appendix ^ 

All the selected forces perform well concerning the 
total energy and radii. They all have comparable in- 
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compressibity if =210-250 MeV and comparable surface 
energy which results from a careful fit to ground-state 
properties. Variations occur for properties which are not 
fixed precisely by ground-state characteristics. The ef- 
fective nucleon mass is 1 for SkT6 and SkP, 0.9 for SkO 
and SkO', around 0.8 for SkM* and Zg-, and even lower, 
around 0.65, for SLy4, Skll, SkI3, and SkM. Isovector 
properties also exhibit large variations. The asymmetry 
energy ranges from very low, 26 MeV for Zo-, to rather 
high, 38 MeV for Skll, with the values for other forces 
being around 30-32 MeV. The appropriate options fo r the 
center-of-mass correction ( |3.7| ) and spin-orbit force ( |3.6| ) 
are found in Table | in Appendix |^. The choice embraces 
^CM=0 as well as 9cm=^, and the various options for 
the spin-orbit force. The only forces in the sample which 
have not yet been published elsewhere are SkO and SkO'. 
They stem from an ongoing exploration of Skyrme forces 
trying to accommodate more observables. In addition to 
SkM which fits ground-state energies, radii, surface thick- 
nesses and the isotope shifts of r.m.s. radii, SkO and 
SkO' also manage to reproduce the jump in the isotopic 
trend of the two-neutron separation energies in the lead 
isotopes, a feature where most Skyrme forces fail. More- 
over, these two forces represent a most recent update of 
the fits along the line of l6^ j60[l now using an up-to-date 

Last but not least, 
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treatment of pairing, see Sec. 
we have here a pair of forces which are fitted in precisely 
the same manner and differ only in the spin-orbit factor 
6\s. This allows for testing specifically the impact of this 
variation. 



B. Treatment of pairing 

In the original publications, various forces were used 
with different pairing recipes. Most of these recipes are 
very schematic (e.g., constant gap or seniority force) and 
fail when proceeding into the regime of exotic nuclei. (See 
discussion in Refs. ||6^,^.) On the other hand, details of 
the actual pairing recipe do not affect the overall quality 
of the forces because these are usually fitted to proper- 
ties of well-bound nuclei. In this work, we compute the 
pairing matrix elements from a local interaction. Among 
several choices, we take the simple S force, which leads 
to the pairing energy functional 



^pair 



E 



^ /dVx^(r)F(r), 



(3.8) 



where the local pair density x(r) reads ||67| , |69|| 

^^UaVa \^a{r)f , (3.9) 



and the function F—1 or F—1 — p{r) / pc gives the volume 
or surface type of pairing correlations, respectively, while 
Pc=0.16fm~^ is the saturation density. 



A further key quantity for pairing is the selection of 
the pairing phase space. Following Refs. |7^,^, to cut 
the space at the high-energy side, we use a Fermi- type 
form factor 



/a 
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l + cxp[(e„ - A- A£;)^]' 



(3.10) 



where the width of the smooth cut-off is linked to the 
energy offset by 
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(3.11) 



as was done in the earlier proposals ||7^,Q. We adopt 
the point of view that pairing is a valence-particle effect 
for which the energy range Ai? should be proportional 
to the average level spacing for the given nucleus and 
nucleon type. We accomplish this by fixing the number 
of pairing-active states: 



: A/'+ 1.65A/'2/^ 



(3.12) 



where Af^Z or N. Details of this "soft" cut-off scheme 
and the reas oning behind the actual choice of the factor 
1.65 in Eq. ( [3.12D can be found in Ref. @. 

Having selected the pairing recipe, one needs to fix the 
strengths Vp and Vn- They have been fitted to empiri- 
cal pairing gaps in a selection of nuclei. Details and the 
resulting pairing force parameters are given in Table || 
in Appendix ^ The surface-pairing strengths for neu- 
trons and protons have been adjusted in the same way as 
the strength parameters of the standard (volume) delta 
pairing. 



C. Description of calculated quantities 

After solving the HF equations in the usual manner, 
we obtain the self-consistent single-particle orbitals from 
which the total energy, as well as several other observ- 
ables, can be calculated, as described in this section. 

A global characteristic of the Skyrme interaction is the 
surface energy coefficient: 



lim 

A-*oo 



\E{A) 






A 


(!) 


A=oo- 



(3.13) 



where E{A) is the energy computed for A particles. Be- 
cause Us is computed from the HF results for large parti- 
cle numbers, it is independent of shell effects, and hence 
it characterizes the surface properties of the bulk energy 
i?ofEq. (1.1). As the limiting process in Eq. (3.13) is ex- 
tremely slow fz^ , it is best to evaluate for semi-infinite 
nuclear matter, and for that we use the semiclassical M. 
Brack code ||]. 

Figure |^ displays the surface energy coefficient for the 
Skyrme parameterizations employed in this work. The 



5 



larger the a^, the greater the surface tension. Conse- 
quently, large values of imply the stronger resistance 
of the system against surface distortions (or, in other 
words, reduced deformability) . As seen in Fig. 0, the 
"stiffest" interactions are SkT6, SLy4, and SkP, and the 
"softest" parameterizations are Z^., SkO, and SkO'. 

A large part of our survey below deals with quadrupole 
deformation potentials. We produce a systematic series 
of deformed mean-field states by adding a quadrupole 
constraint Q cx r^l2o/[l + /(^)] to the HF field, where 
the function /(r) suppresses Q at large distances (see 
Ref. 0). The calculated deformed shapes are charac- 
terized by means of the dimcnsionless quadrupole defor- 
mation: 



The rotational moment of inertia is determined from 



5 A{r 



(3.14) 



The total energy of the system iJtot as a function of (3 
represents a zero-order approximation to the potential 
energy curve for /3-vibrations, i.e.. 



v{(3) = Etom = {^p\m0). 



(3.15) 



However, before one can use V{(3) in calculations with the 
collective Hamiltonian, dynamical corrections have to be 
added. The reason is that the underlying states |$^) 
have a finite uncertainty in the collective deformation, 
i.e., ^0. As a consequence, the potential V{f3) con- 
tains contributions from /3-fluctuations in \^i3), and these 
contributions need to be subtracted first before adding 
the energies associated with the true physical zero-point 
fiuctuations in (3. The theoretical evaluation of these cor- 
rection terms can be done in the framework of the Gen- 
erator Coordinate Method at the level of the Gaussian 
Overlap Approximation (GOA), as has been discussed in 
several publications. (See Ref. [|2j for a review.) The 
collective parameters in the present (axially symmetric) 
case are the quadrupole deformation [3 and the two ro- 
tational angles •& and tf. The volume element in these 
coordinates is not Cartesian and thus one has to em- 
ploy the GOA in a topologically invariant fashion. For a 
detailed discussion of the general case, see Ref. |74 . Sim- 



pler formulae used in this work are taken from Ref. [[ll[ 
namely, we define the corrected deformation energy as 



V{(3) = V{f3) - {EzPE.p + £^ZPE,rot) , 



(3.16) 



where the rotational, i?zPE,rot, and vibrational, i?zPE,/3, 
zero-point energy corrections read 



ZPE.rot 



ZPE,/3 



2 6 rot 

2{dpHdii)-(dlH)-{Ht 



/3/ 



3-.g 



-£^ZPE,io 





= {[L,[H,L]]), 



(3.19) 



while the switching factor g ^(/^)/2^, which originates 

from the topologically invariant extension of GOA pT] , 
is defined as 



/o da;?7(x2-l)e''("'-i) 



(3.20) 



In Eqs. (p^)-(p^),the average values () are taken with 
respect to the /3-dependent HF states |<&/3). The defini- 
tion of the collective mass parameters recurs, in princi- 
ple, to the full Hamiltonian H. However, for the present 
exploratory purposes, we employ the Inglis cranking ap- 
proximation which is obtained from the above expres- 
sions by letting H ^ ho, with ho being the mean-field 
Hamiltonian. In the following, the results of calculations 
of the potential energy surfaces (PES) always pertain to 
the tot al ene rgies corrected for the zero-point motion, as 
in Eq. (|T6|). 



D. Discussion of Potential Energy Surfaces 

F or the set of Skyrme parameterizations described in 
Sec. |IIIA|, the PESs have been calculated for 26,28,30,32]s^g^ 



30,32 34 jyjg 38,40,42,44g 80,82,842^. g^j^^J 92,94,96,98,1002j- 

functions of quadrupole deformation /3. These results are 
discussed below. 



1. Deformation in the N~20 region 

The results of calculations for 30,32,34]yjg g^^.^ shown in 
Fig. ^. For most Skyrme parameterizations used, the 
pattern is fairly similar. Namely, the nucleus '^"Mg is 
predicted to be merely deformation-soft, while the oc- 
cupation of the fi/2 neutron shell in ■^^Mg gives rise 
to a very deformed intrinsic shape with P ranging from 
0.3 to 0.4. The nucleus "^^Mg appears to be a transi- 
tional system with coexisting spherical and prolate min- 
ima. For Skyrme parameterizations SkM* and Z^, the 
prolate minimum is calculated to be practically degen- 
erate with the spherical one. For the remaining forces, 
the prolate structure (sometimes corresponding to a local 
minimum, sometimes forming a shoulder in the PES) lies 
from 2MeV to 4MeV above the spherical ground state, 
depending on the choice of the Skyrme parameterization. 

A similar pattern is observed in Fig. B| for the neutron- 
rich Ne isotopes. Here, the nuclei ^^'^°Ne are predicted 
to be very soft, strongly anharmonic, while "^^Ne is well 
deformed in all cases. The semi-magic "^^Ne is predicted 
to be spherical. However, as in the case of "^^Mg, a low- 
lying secondary prolate minimum develops in the SkM* 
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and Z„ models. By comparing Figs. || and || one notices 
that for all the forces used, the deformed configuration 
in ^°Ne lies ~1 MeV higher in energy than that in ^^Mg. 
That is, the shape mixing phenomenon is expected to be 
much stronger in 32 Mg than in ^"Nc. 

Of course, in the case of low-lying coexisting states, the 
energy difference between spherical and deformed min- 
ima depends strongly on the details of the calculations. 
In particular, variations in the treatment of pairing cor- 
relations are expected to play a role in light nuclei such 
as 32 Mg. To illustrate this point, we performed two addi- 
tional sets of calculations for ^^Mg using different pairing 
recipes. Figure I shows the PESs for 30,32,34 j^jg obtained 
by taking (i) volume pairing as in Fig. B, (ii ) the sur- 
face pairing interaction as defined in Eg. (|3.8| ), and (iii) 
neglecting pairing (i.e., pure HF). 

As expected, the prolate minimum is well developed 
in most unpaired calculations, and its energy is signifi- 
cantly lowered as compared to the calculations with pair- 
ing. (For the forces SkM*, and Skll the prolate un- 
paired minimum becomes the ground state.) The oppo- 
site holds for the surface-pairing variant: the correspond- 
ing PESs seem softer in the direction of /3. The sensitivity 
of the calculated excitation energy of the intruder state 
in 32Mg on the pairing recipe indicates that the detailed 
description would require (i) a realistic pairing interac- 
tion that could be applied in mean-field calculations for 
light nuclei, and (ii) the proper treatment of particle- 
number fluctuations. Other uncertainties in determining 
the relative energies of coexisting states are discussed in 
Sec. [HE below. 



There are many factors that can influence the energy 
difference between coexisting states. Probably the most 
important one is the single-particle shell structure. Posi- 
tions of individual shells are strongly affected by changes 
in Skyrme parameters, in particular those defining the 
spin-orbit term. 

The spherical neutron shell structure for ^^jyjg pre- 
dicted by various Skyrme parameterizations is shown in 
Fig. H Of particular interest is the size of the A^=20 
magic gap which is measured by the distance between 
the /7/2 and (^3/2 shells: 



Ae2o = e(/7/2) - 6(^3/2). 



(3.21) 



The variations of Ae2o are nicely correlated with the be- 
havior of the height of the prolate minima AiJpioi in 
^^Mg, shown in Fig. ^ for two variants of calculations: 
with and without pairing (the latter to single out the pure 
effect of the particle- hole channel). Indeed, the large val- 
ues of AiJproi in Skis, SkM, and SkO can be correlated 
with large values of Ae2o. Likewise, small shell gaps in 
SkM* and Z^- are consistent with AE'pi.oi^O obtained in 
these models. However, there are exceptions to this rule. 
For instance, the value of Ae2o is rather low in SkT6 but 
the prolate minimum is calculated to be at ^3 MeV. 

In order to better understand some of the deviations 
between the pattern of Ae2o and AiSproi, it is instructive 



to return to Fig. 0. The "stiffest" interactions are SkT6, 
SLy4, and SkP, and - indeed - for all these forces, spher- 
ical ground states are predicted. The "softest" param- 
eterizations are Z^-, SkO, and SkO', but the large value 
of Ae2o in SkO and SkO' gives rise to spherical ground 
states. 

The summary of single-neutron ^3/2 and /7/2 energies 
for the A^=20 isotones calculated with several Skyrme 
forces is shown in Fig. ||. As expected, the absolute 
binding energy of these shells decreases rapidly when ap- 
proaching the drip-line nucleus ^^O. For all the interac- 
tions considered, however, Ae2o varies very slowly with 
Z. 



2. Deformation in the N^28 region 

The results of calculations for 38,40,42,44g shown in 
Fig. I indicate that the A^=28 shell gap is broken around 
^^S. Indeed, most interactions used predict a deformed 
ground state for ^*S. It is worth noting that the two pa- 
rameterizations that yield strongest deformation effects 
in ^^Mg, namely SkM* and Z^y, do not produce deformed 
minima in '^'*S but rather /3-unstable PESs. This indi- 
cates that the deviations between results should be linked 
to the details of the underlying shell structure which 
looks, of course, different for the different shell closures. 

Figure || shows the single-neutron structure in "'^S to- 
gether with the calculated energies of the prolate, Ai?p_s, 
and oblate, AE'ojS, minima (with respect to the spherical 
configuration). The position of the deformed minimum 
is greatly influenced by the size of the N=28 gap 



Ae28 = e(p3/2) - e(/7/2). 



(3.22) 



For most interactions considered, Ae28 is small - of the 
order of 2-3 MeV. Consequently, in most cases, the de- 
formation energies follow the pattern of Os. 



3. The N^Z^40 region 

The interplay between spherical and deformed subshell 
closures at N or Z=40 is illustrated in Fig. 0. Although 
coexisting spherical and prolate minima in *''Zr are pre- 
dicted for all the Skyrme parameterizations used, their 
relative position does depend strongly on the interaction. 
The interactions SkM*, Z„, Skll, SkM, and SkO' predict 
a strongly deformed ground state for ^"Zr, in agreement 
with experiment. Other forces, most notably SkP and 
SkT6, yield a spherical ground state. 

The spherical shell structure in ^°Zr is displayed in 
Fig. Since for this nucleus Z—N, the proton and 
neutron single-particle energies are very similar. (The 
influence of Coulomb interaction on shell structure in 
this medium-mass system is weak.) As in the case of 
32 Mg, there is a clear correlation between the size of the 
N=Z=AO subshell closure, 
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Ae4o = e{gg/2) - e(pi/2), 



(3.23) 



the deformation energy, and the surface-energy coeffi- 
cients. For all Skyrme parameterizations which predict a 
spherical ground state in ^°Zr, either Ae4o is large (like 
in SkI3) or is large (like in SkP), or both. 

The PES and corresponding shell structure of ^°Zr pro- 
vide a particularly clear example of how variations in the 
treatment of the spin-orbit force can have a large im- 
pact on the results. Compare the "twin" parameteriza- 
tions SkO and SkO' which differ just by the switch in 
the spin-orbit functional (3.6). The different spin-orbit 
force produces a different splitting of the Ig levels, sub- 
sequently a different shell gap at the Fermi surface (see 
Fig. and finally a different excitation energy of the 
prolate minimum (see Fig. ^ and the PES in Fig. 



4. The N^56, Z^40 region 

In this region of shape coexistence, the best agreement 
with the observed experimental trend is given by SkM*, 
SLy4, and Skll (see Fig. Namely, ^^Zr is pre- 
dicted to be spherical, -'^'^^Zr very well deformed, and ^^Zr 
spherical, with a low-lying deformed intruder state. The 
worst agreement with the data is obtained in the SkP 
model in which all isotopes considered have spherical 
ground states, and in the SkI4 model which predicts a 
strongly deformed ground state for 94,96,982j-^ 

Again, the general pattern of deformation energies can 
be explained in terms of the calculated gap sizes: the 
Ae4o proton gap and the A'^=56 gap 



Ae56 = e(si/2) - 6(^5/2). 



(3.24) 



For instance, for the interaction SkM the proton Ae4o 
and the neutron Aese are rather small and this yields a 
deformed ground state in ^^Zr. The opposite holds for 
SLy4, which, in addition, has a large value of as- Hence, 
it predicts spherical ^^Zr. 



E. Zero-point fluctuations 



reduced, and this causes -EzPE,rot to drop. The differ- 
ence of £'zpE,rot between spherical and deformed minima 
is around 4 MeV, i.e., this is a significant correction to the 
total energy. As discussed in Ref. however, the rota- 
tional zero-point energy should be su pplem ented by the 
vibrational counterpart Ezpe./s, Eq. ( 3.18| ). This quan- 
tity shows an opposite behavior: it is strongly peaked 
around the spherical shape and reaches the value of 
~1 MeV at large deformations. The large peak at zero de- 
formation compensates for the correspondingly large dip 
in the rotational ZPE such that, altogether, a smooth to- 
tal ZPE emerges whose main variation is the global trend 
to grow with deformation. The irregularities (kinks) in 
EzPE,[3, seen in Figs. p^[]T6| , are caused by the unphys- 
ical collapse of the BCS pairing in certain regions of f3, 
which, in turn, produces enormous spikes in the collec- 
tive quadrupole mass. Clearly, it is necessary to improve 
the description of zero-point fluctuations by (i) taking 
into account the particle-number fluctuations, and (ii) by 
going beyond the Inglis cranking approximation. Based 
on the present results, however, one can conclude that 
the zero-point correction should be rather small for '^^Mg 
and ^''S, and that it favors the deformed state by about 
2 MeV for ^OZr and about IMeV for ^^Zr. The effect of 
shape fluctuations becomes more important at large de- 
formations due to the steady increase of -EzPE.rot- Con- 
sequently, when studying superdeformations, fission bar- 
riers, fission valleys, etc., zero-point corrections should 
be taken into account. 



IV. SHELL-MODEL ANALYSIS 

The mean-field analysis presented in the previous sec- 
tion is supplemented by shell-model calculations for the 
neutron-rich nuclei around ^^Mg using the Shell Model 
Monte Carlo (SMMC) technique |^,^. In contrast to 
the mean-field approach, shell-model calculations prop- 
erly treat configuration mixing and dynamical fluctua- 
tions. On the other hand, the rather small conflguration 
space employed (here, two oscillator shells) in compari- 
son to the mean fleld can lead to an improper description 
of certain states. 



The role of fluctuations beyond the mean fleld is illus- 
trated in Figs. p^[]T6| which show the effect of rotational, 
/3-vibrational, and center-of-mass corrections. The calcu- 
lations were performed with the SkM parameterization; 
a very similar result (not shown here) was obtained with 
the SkM* force. 

The center-of-mass correction, Eq. { \i.7\) , depends very 
weakly on deformation; hence its contribution to the de- 
formation energy can be safely neglected. The rotational 
zero-point energy, Eq. (3.17), is zero at the spherical 
shape and increases steadily with deformation. The ad- 
ditional fluctuations of £'zPE,rot with P are mainly due 
to the changes in the pairing fleld: the moment of inertia 
Grot, Eq. (3.19), increases when pairing correlations are 



A. Shell-Model Monte Carlo Method 

The SMMC method offers an alternative way to calcu- 
late nuclear structure properties, and is complementary 
to direct diagonalization. SMMC cannot, nor is it de- 
signed to, flnd every energy eigenvalue of the Hamilto- 
nian. Instead, it is designed to give thermal or ground- 
state expectation values for various one- and two-body 
operators. Indeed, for larger nuclei, SMMC is presently 
the only way to obtain information on properties of the 
system from a shell-model perspective. 
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The partition function of the imaginary-time many- 
body propagator, U = exp{—f3H), is used to calculate 
the expectation values of any observable 



TrUn 
TrU 



where 



H = Hi+H2 



(4.1) 



(4.2) 



is the shell-model Hamiltonian containing one-body and 
two-body terms, and T—\/ [3 is the temperature of the 
system. The two-body term, H2, is linearized through 
the Hubbard-Stratonovich transformation, which intro- 
duces auxiliary fields over which one must integrate to 
obtain physical answers. Since H contains many terms 
that do not commute, one must discretize (3 = NtA(3. 
The method can be summarized as 



Z = TrU = Trexp(-/3ff) ^ Tr exp(-A/3ff) 



P[cr]G((T)Tr Yl exp A/3/i(cr„) 



n=l 



(4.3) 



where cr„ are the auxiliary fields. (There is one cr-field for 
each two-body matrix-element in H2 when the two-body 
terms are recast in quadratic form.) I?[cr] is the measure 
of the integrand, G(a) is a Gaussian in a, and ft, is a 
one-body Hamiltonian. Thus, the shell-model problem is 
transformed from the diagonalization of a large matrix 
to one of large dimensional quadrature. Dimensions of 
the integral can reach up to 5-10* for the sdpf systems, 
and it is thus natural to use Metropolis random walk 
methods to sample the space. Such integration can most 
efficiently be performed on massively parallel computers. 
Further details are discussed in Ref. | |76| |. 

The SMMC method is not free of extrapolation when 
realistic Hamiltonians are used. The sign problem for re- 
alistic interactions was solved by breaking the two-body 
interaction into "good" (without a sign problem) and 
"bad" (with a sign problem) parts: H = Hgood + Htad- 
The part Hbad is multiplied by a parameter, g, with val- 
ues typically lying in the range — 1 < <? < 0. The Hamil- 
tonian H{g) = f{g)Hgood + gHbad has no sign problem 
for g in this range. The function f{g) is used to help in 
extrapolations. It is constructed such that f{g = 
and takes the form [1 — (1 — g)/x]i with x — 
The SMMC observables are evaluated for a number of 
different negative g-values, and the true observables are 
obtained by extrapolation to g=l. A prescription has 
been used to remove center-of-mass contaminations in- 
herent in the wave functions when multi-?iw spaces are 
used . In each calculation presented here, we took 6 
values of g, and 4096 independent Monte Carlo samples 
per 5 value. 



B. The Effective Shell-Model Interaction 

In this work we wish to compare two shell-model in- 
teractions that could prove useful for the sdpf region. 
The first interaction was derived using microscopic tech- 
niques j?^ , while the second is a more piece- wise inter- 
action similar to those used in highly truncated standard 
shell-model calculations for nuclei near N=20. 

Our first interaction, dubbed sdpf, is described in de- 
tail in Ref. |^ . In order to obtain a microscopic effective 
interaction, one begins with a free nucleon-nucleon inter- 
action V which is appropriate for nuclear physics at low 
and intermediate energies. The choice made in Ref. |79| ] 
was to work with the charge-dependent version of the 
Bonn potential models as found in Ref. |^^. Standard 
perturbation techniques were then employed to obtain an 
effective interaction in the full sdpf model space. The in- 
teraction was then modified in the monopole terms using 
techniques developed by Zuker and co-workers ]8l| , |8^ . 

The second shell-model interaction employed in this 
work, dubbed sdpf, results from a more standard, yet 
less rigorous, approach to the problem. Numerous shell- 
model studies have been carried out in truncated model 
spaces for neutron-rich nuclei near N—20 ]2^-]2^ and 
7V=28 [||,|9|j2|]. Several sdpf shell effective interac- 
tions were used in these studies; many of these inter- 
actions are quite similar in a number of respects. All of 
them use the Wildenthal USD interaction |^3| in the sd 
part of the Hilbert space. All also use some 'enhanced' 
version of the original Kuo-Brown p/-shell G-matrix in- 
teraction |84j to describe nucleons in that shell. The 
cross-shell interaction is handled in one of two different 
ways: matrix elements are generated via a G-matrix or 
via the Millener-Kurath potential |85|. As is common 
in this type of calculation, selected two-body matrix ele- 
ments and single-particle energies have been further ad- 
justed to obtain agreement with experiment. Here, we 
use the following prescription: we incorporate the USD 
interaction for the sd-shell |^], and the FPKB3 inter- 
action as found in Ref. |Q. We also used the standard 
Millener-Kurath |^ prescription for the cross-shell ma- 
trix elements. However, our first investigations found 
that the scattering of particles from the sc?-shell to the 
p/-shell was too strong. Therefore, we reduced the cross- 
shell monopole matrix elements by 1.4 MeV. The single- 
particle energies were adjusted to fit ^^Ca single-particle 
energies 

The sdpf interaction describes satisfactorily the 
ground-state masses in the sd-pf region. The differ- 
ence between theory and experiement in the binding en- 
ergies for the 10 nuclei studied in Ref. ||7^ is approx- 
imately ±1.5 MeV with a statistical error of 0.75 MeV. 
B{E2) values were well described across the sd-pf region 
using standard effective charges (ep=1.5 and e„=0.5). 
Occupation probabilities for the /7/2 shell were in fair 
agreement with highly truncated interaction scenarios. 
The sdpf interaction cannot describe the B{E2) values 
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across the sd-pf region unless one invokes two sets of 
effective charges (ep=1.5, e„=0.3 in the A<AO region, 
and ep=1.2, e„=0.1 in the A>40 region). Furthermore, 
binding energies were not weU reproduced in the sdpf 
interaction, although the excitation spectrum for a light 
nucleus (e.g., ^^Mg) was of the same quality as that of 
the sdpf interaction. The occupation of the full p/-shell 
in the neutron-rich nuclei is similar in both the sdpf and 
sdpf interactions by construction, although more par- 
ticles occupy levels other than /7/2 in the sdpf case. 
B{E2) values and occupations numbers of three nuclei 
were used in the fitting procedure of sdpf: ^^Ar, ^^Mg, 
and *'*Ti. Thus, it is not surprising that the behavior of 
the two interactions is similiar around '^^Mg, while dif- 
ferences occur for other nuclei (see discussion below). 

It should be clear that we prefer the sdpf interaction 
as it is based more on a theoretical derivation across the 
entire shell-model space in which the calculations were 
performed. However, we believe it is worthwhile to inves- 
tigate the differences between this interaction and those 
obtained in a more phenomenological way, such as sdpf . 
We also note that interactions derived in a similar fashion 
to sdpf have served very useful purposes when calcula- 
tions using them are performed in truncated spaces (e.g., 
as those by Retamosa et al. |^^). However, they are less 
able to reproduce experimental data in full-space calcu- 
lations such as those performed here. 



C. Results of Shell-Model Calculations 

The SMMC calculations were performed for a number 
of even-even nuclei from the neutron-rich iV=20 region. 
In order to relate the SMMC results to the schematic 
shell-model scheme based on the broken-pair approach 
0,|8|, we show in Fig. |l^ the mean value of QpQn, 
related to the pr oton - neutron quadrupole interaction en- 
ergy Eq of Eq. (l.S), and the mean value of ^J^Aqi, re- 
lated to the pairing energy impair in the J=0, T~l channel 



(Agi is the J=0, T=l pair operator Q). The calcula- 
tions were performed for the neutron-rich Ne and Mg iso- 
topes. The corresponding orbital occupation coefhcients. 



2j„ + l' 



(4.4) 



where Nj^ is the average number of particles in the shell 
ja , are displayed in Fig. [l8| . 

For the sdpf interaction, the result is consistent with 
the trend predicted by the schematic model. Namely, the 
expectation value of QpQn increases at A^=20 and 22, 
reflecting the increased occupation of the fj /2 shell. For 
the sdpf interaction, however, the pattern is markedly 
different. In particular, (QpQn) varies very little with 
N, especially for the Mg isotopes. Although the sdpf 
interaction predicts larger occupations of the /7/2 shell, 
the value of Eq seems to be significantly greater in the 



sdpf c ase. W e shall come back to this apparent paradox 
in Sec. |IVD| . 

Both interactions yield fairly constant (A^-^Aqi) for the 
protons (the proton pairing energy does not change with 
neutron number) and an almost linear increase with N 
for the neutrons (this behavior is indicative of a weak 
neutron pairing). In order to understand an extremely 
weak dependence of neutron (A^-, Aqi) predicted in the 
sdpf calculations, we show in Fig. the J=0, T=l ma- 
trix elements, {jaja01\V\jf}jpOl) , of sdpf and sdpf. It 
is seen that, in general, the pairing interaction within the 
sd and fp shells is weaker for sdpf, and the opposite is 
true for the cross-shell pair scattering. Moreover, except 
for the ^5/2 shell, the diagonal pairing matrix elements 
{a=(3) of sdpf are either close to zero or positive (i.e., the 
pairing interaction in these states is actually repulsive!). 



D. Mean- field analysis of shell- model results 

The shell-model Hamiltonian (4.2) can be written as 
H = eaaitta + ^ X! ^"/37'5 aaajjO^a^, (4.5) 

where the single-particle indices (indicated by Greek 
letters) denote the single-particle quantum numbers 
(n,l, j^m^T = tz), Ca are the single-particle shell-model 
energies, and Vap-^s are the (antisymmetrized) two-body 
matrix elements of the two-body interaction. 

In order to translate shell-model results to the language 
of mean-field theory, we carried out the HFB calculations 
using the shell-model Hamiltonian (4.5). In the following, 
this variant of calculations will be referred to as HFB- 
SM. In the calculations we impose spherical symmetry 
and disregard neutron-proton pairing. The details of the 
HFB-SM derivations are given in Appendix 

The canonical HF B si ngle-particle d^/2 and f^ /2 neu- 
tron energy levels ( Blj) calculated in the sdpf vari- 
ig. |20| for the Ne and Si isotopes. 



ant are shown in Fig. 
Based on this result, two interesting conclusions can be 
drawn. First, the isotonic dependence of single-particle 
levels is very weak. Consequently, the size of the 7V=20 
gap varies little with N (this conclusion also holds for 
the sdpf interaction). Second, the single-particle ener- 
gies strongly depend on Z . This effect has been noticed 
in Ref. p^ , and was discussed therein in terms of the 
monopole neutron-proton interaction, that is, the shift 
in the spherical single-particle neutron energies due to 
protons. It is seen that this monopole effect gives rise 
to the reduction of the iV=20 gap when decreasing Z. 
Indeed, as shown in Fig. |2^, the size of the iV=20 neu- 
tron gap calculated with the sdpf interaction decreases 
from -lOMeV in to ~2MeV in ^sq. It is impor- 
tant to emphasize that this monopole effect predicted in 
HFB-SM, important for the excitation energy of the de- 
formed intruder configuration p8,p[, is not a threshold 
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phenomenon due to the weak binding; the reduction of 
the magic gap comes solely from the shell-model interac- 
tion. 

Figure also shows the value of Ae2o predicted with 
the sdpf interaction. Here, the dependence of the gap 
on the neutron number is very weak. To understand 
the difference between predictions of the two interacti ons , 
Fig. ^ shows the matrix elements V{a,l3) of Eq. (BE) 
for sdpf and sdpf . These particle-hole matrix elements 
define the self-consistent mean-field, hence the canonical 
single-particle energies. Since the single-particle shell- 
model energies do not vary with particle number, the 
variations of Ca with TV and Z are solely due to changes 
of the self-consistent mean-field. In addition, since the 
neutron-neutron contributions to Ae2o do not depend on 
Z, the variation of the iV— 20 gap with proton number 
can be traced ba ck t o the proton-neutron interaction. 
According to Eq. (B4), the main contribution to the Z- 
dependent part of Ae2o comes from the proton-neutron 
terms: 



Ernrr — E[ 



SM 



HFB- 



(4.7) 



(4.6) 



For nuclei discussed in Fig. ^the occupied proton shells 
are ^5/2 a nd S i/2, and, precisely for these orbitals, the 
difference (4.6) is close to zero for sdpf and it is about 
1 MeV for sdpf . That is, it is very close to what is seen 
in Fig. ^ (Ae2o changes by ~1 MeV/proton). One can 
thus conclude that the monopole effect of Refs. IH,^ is 
very weak for the spdf interaction. 

Coming back to the prediction of the sdpf interaction 
concernin g the unexpected behavior of {QpQn) versus N 
(see Sec. [V C ), it is instructive to inspect the particle- 
hole matrix elements of Fig. The proton-neutron ma- 
trix elements of sdpf are negative (i.e., the particle-hole 
interaction is attractive in this channel), and they are 
significantly larger in magnitude than the like-particle 
matrix elements (the latter ones are usually attractive 
or close to zero). This result does not come as a sur- 
prise; it is generally believed that the proton-neutron 
component of the particle-hole interaction is dominant 
|p|, ^7|j5^ . For sdpf, however, the situation is differ- 
ent: the proton-neutron interaction is, generally, much 
stronger (especially in the pf shell and for the cross- 
shell matrix elements), but the particle-like matrix ele- 
ments are all positive. Therefore, the structures predicted 
in sdpf result from a subtle balance between strongly 
attractive proton-neutron particle-hole interaction and 
repulsive (and weaker, see Fig. p2|) proton-proton and 
neutron-neutron particle-hole forces. This is reflected in 
the SMMC results shown in Tabic ||. In sdpf, the values 
of {QpQn) and (Qn) steadily increase when crossing the 
N—20 gap, consistent with the increasing f'j/2 occupa- 
tions. This is not the case for sdpf where the quadrupole 
collectivity decreases in '^''Mg in spite of the fact that the 
/7/2 occupations are larger and the Ae2o gap is smaller 
than in the sdpf model. 

Figure p3| shows the SM correlation energy, i.e., the ex- 
cess of binding energy above the spherical HFB estimate 



For sdpf the behavior of i?corr does not follow the pat- 
tern of increased quadrupole collectivity when crossing 
the N=20 gap. Actually, i?corr decreases. This result is 
consistent with our mean-field results which predict the 
coexistence of spherical and deformed shapes in '^^'•^'*Mg, 
giving rise to the quadrupole-softness or shape mixing. 
Again, the behavior of -Ecorr in sdpf is different. There 
is very little change in the SM correlation energy for the 
Mg isotopes; its rather large value reflects the increased 
correlations due to the significant occupation of the f-^ ji 
shell (spherical HFB-SM calculations predict no pf neu- 
trons in ■^^Mg). For both forces, the correlation energy 
in ■^^Mg is greater than in ■^"Ne. This result corrobo- 
rates our HF prediction that "^^Ne is more spherical (i.e., 
coexistence effects are weaker) . 

To see the sensitivity of the SMMC predictions for 
■^^Mg to the size of the splitting between the undisturbed 
single-particle energies Cq, we changed the splitting by 
±0.5 MeV around the standard value. Surprisingly, such 
a variation changes the B{E2) (or {Q'^)) value and the 
shell-model occupation coefficients very little. The corre- 
lation energy changes from -16.9 MeV (standard N=20 
sphtting) to -16.6 (shell gap decreased by 0.5 MeV) and - 
12.4 MeV (shell gap increased by 0.5 MeV). Hence, again, 
in the sdpf interaction the correlation energy is not ob- 
viously related to the quadrupole collectivity. 



V. CONCLUSIONS 

In this paper, we have studied the phenomenon of 
shape coexistence in semi-magic Mg-, S-, and Zr-isotopes 
employing two complementary theoretical approaches, a 
self-consistent mean-field model (Skyrme-Hartree-Fock) 
and shell-model calculations which account for all cor- 
relations in a restricted space. The main conclusions of 
this study can be summarized as follows. 

The variety of Skyrme-HF predictions has been ex- 
plored by comparing all results for a set of 10 typical 
effective Skyrme forces. For mean-field models, shape 
coexistence can be quantified in terms of the relative en- 
ergies of coexisting local minima. All selected Skyrme 
forces agree in producing the same isotopic trends in 
these key features of shape coexistence, but the actual 
preference for a spherical or deformed ground state varies 
from force to force. We have tried to relate the results 
to other important features of the nucleus and find that 
the main factor that determines the excitation energy of 
the deformed intruder state in the HF calculations is the 
single-particle shell structure (in particular, the sizes of 
the spherical magic gaps and subshells). Another im- 
portant quantity that defines the nuclear deformability 
is the surface energy coefficient as. Skyrme interactions 
with large values of ag (SkT6, SLy4, SkP) favor spherical 
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configurations as compared to other forces (provided that 
the corresponding shell effects are similar) . On the other 
hand, forces with low values of as {Z^, SkO, SkO') give 
rise to softer PES and low-lying intruder configurations. 

The single-particle structure can be strongly affected 
by small variations in the definition of the energy func- 
tional. In this context, a good example is the treatment 
of the spin-orbit term by various parameterizations with 
respect to the inclusion of the contribution. For this 
purpose we had a twin pair of forces (SkO and SkO') 
in the sample which differs just by this feature. It was 
found that this modification can have a large impact on 
shape coexistence in some cases (here the most dramatic 



is 



80 



Zr). 



The proper treatment of pairing and zero-point cor- 
relations is crucial if one aims at detailed predictions of 
shape coexistence. For instance, according to our es- 
timates, the zero-point rotational-vibrational correction 
should be around 2MeV in '^"Zr, around 1 MeV in ^'^Zr, 
and is expected to increase systematically with deforma- 
tion. 

For the Skyrmc interactions considered, the size of the 
N=20 gap varies very slowly with Z, and, except for 
SkT6 and SkP, Ae2o is quenched when approaching ^^0 
(see Figs. ^ and |4|). This result agrees with the sdpf 
HFB-SM calculations. On the other hand, the size of 
the A^=20 neutron gap calculated with the sdpf inter- 
action decreases rapidly with Z. This strong monopolc 
effect can be traced back to differences between certain 
proton-neutron matrix elements of the shell-model inter- 
action |^,|| . It is important to emphasize that this effect 
has its roots in the properties of the shell-model Hamilto- 
nian and should not be confused j27| with the threshold 
phenomena due to weak binding and the closeness of the 
particle continuum. Also, for a given isotopic chain, the 
N dependence of the A^=20 gap has been found very 
weak for both shell-model interactions. This contradicts 
recent conclusions of Rcf. ||2^ which predict the sharp 
minimum of Ae2o at iV=20. It should also be noted that 
the size of the single-particle gap does not always corre- 
spond to the shell-gap parameter related to a difference 
between two-neutron separation energies: 



d2n = S2niN)~S2niN + 2). 



(5.1) 



Indeed, as seen in Fig. based on the spherical Skyrme- 
HFB calculations, while Ae2o changes very weakly with 
iV, d2n experiences a dramatic drop when approach- 
ing Z=8. This indicates strong effects related to self- 
consistency in light drip-line nuclei. 

The nucleus '^^Mg has been found to be a classic ex- 
ample of shape coexistence; the spherical and deformed 
configurations are close in energy and shape mixing is 
expected. This prediction is consistent with the recent 
measurement from GANIL [Q according to which the 
E4+ / E2+ ratio in ^^Mg falls well below the rotational 
limit. A similar mixing effect is predicted to occur also 
in ^°Ne but is much weaker. For most Skyrme pa- 
rameterizations used, the A^=28 gap is predicted to be 



rather small. This gives rise to strong deformation ef- 
fects around ^''S. The strong coexistence effects are also 
predicted for ^°Zr and ^*Zr. 

Both families of models applied in this work, i.e., self- 
consistent mean-field models and the shell model, should 
be viewed as effective theories. That is, their predic- 
tive power crucially depends on the effective interaction 
assumed. Since we do not know the "true" energy func- 
tional (though we know that it exists ]88| , |89| ), and we are 
still unable to derive "exactly" the effective shell-model 
interaction and the effective shell-model operators, we 
are bound to try different parameterizations. From this 
point of view, nuclear coexistence is a very challenging 
battleground. Although the global picture is understood, 
the structural details strongly depend on the actual phe- 
nomenology used and approximations involved. 
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APPENDIX A: THE SKYRME 
PARAMETRIZATIONS 

For completeness, we provide the parameters for the 
sample of ten representative Skyrme forces used in this 
study. The parameters hi and h[ used in the definitions of 
Sec. [II A are chosen to give the most compact formula- 



tion of the energy functional, the corresponding mean- 
field Hamiltonian, and residual interaction. They are 
related to the standard Skyrme parameters ti and Xi 
ill]] by: 



ho 
hi 

hi 

b'o 
b'l 
b'2 
b'. 



= to{l + \xo), 



tx{l + \xi) + t2[l + \X2)\ , 
iti{l + \xi) - t2{l + \X2)\ , 

^3(l + 5a;3), 



(Al) 



to{\ + xo), 



f 



3tl(i 



ri3(,2 



t2{\+X2)\ , 
+ t2(| +2:2)] 



•T3J 



Tabl e |^ displays the parameters of the Skyrme func- 
tional (3.3) give n in the form recoupled to the ti, Xi ac- 
cording to Eq. ( |Al[ ) (most of the existing codes use this 
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form of input). All conventional Skyrme forces used sim- 
pler pairing recipes. The pairing strengths Vj, a nd Vn 
for the present pairing treatment (see Sec. [II B) have 
been adjusted anew to the neutron gaps in iT2TTa3,i24gj^ 
(using the values 1.41, 1.39, and 1.31 MeV respectively) 
and proton gaps in ^^^Xe and ^^''Sm (using 0.98 and 1.25 
MeV). The forces SkO and SkO' contained these gaps in 
the pool of data throughout the fit. 



APPENDIX B: THE HFB APPROXIMATION TO 
THE NUCLEAR SHELL MODEL 

The antisymmetrized t wo- b ody matrix clement of the 
shell-model Hamiltonian (^^) can be written as: 

Vq/37,5 = {jamaTa, jl3mi3Ti3\H2\j^m^T^, jsmsTs) 

= X1("1)^°'^'"'^^^"^''~^'(2J+ 1)(2T + 1) 



JT 



rria mp —M J \ ms —M 



X {.Uj0JT\H2\j^jsJT)^{l + d„p)(l + d-,s), (Bl) 
where the condition 

{jo.JpJT\H2\j^jsJT) = i-iy+T-^^-^^ 

X {jc.jpJT\H2\j5jyJT) (B2) 

guarantees the antisymmetrization of matrix elements. 

Because of the condition of sphericity, and the fact 
that in the shell-model space considered each spherical 
shell has a unique value of {l,j), the HFB procedure is 
particularly simple. Namely, the quasiparticle canonical 
states are given by a BCS transformation 



(B3) 



The amplitudes (wq,Wq) define the self-consistent mean 
field: 



(B4) 
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the self-consistent pairing gaps 

/3 



X v'2i/3 + l(jaJa01|il2|j/3j/301), 

and the total HFB energy: 



(B5) 



(B6) 



In deriving Eq. (B5) we employed the phase convention 



of Condon-Shortley for time reversal: 

f\nljm) = (-l)'+^'+"|nZj - m). 



(B7) 



The particle-hole matrix element Vafjafj in Eq- (B4) can 
be written as 



where 



1 + 6a(3 



(B8) 



(2j, + l)(2j^ + l) 
and the (2 J -I- l)-averaged matrix elements are 



(B9) 



^l7^ = E(2'^ + l)(jW/3^1|^2|jaJ/3Jl), (BlOa) 



^l7^ = T.^^J + l)(j„j/3 J0|7?2baj/3 JO). (BlOb) 
J 

The neutron and proton Fermi levels, Xt^ , are determined 
from the particle number equations 



(Bll) 



where N1/2 and N_i/2 are the numbers of valence neu- 
trons and protons, respectively. The HFB equations are 
reduced to a set of coupled equations for occupation am- 
plitudes: 



1 



A* 



\/(ea - At J2 -t- A2 



where 



-■a ^ 



(B12) 



(B13) 



are c anon ical single-particle energies. Equations (Bll) 
and (|B12| ) have been solved iteratively. 
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TA BLE I. Parameters of the Skyrme forces used in this study given in terms of the functional as specified in Sec. Ill A and 
[II E|. The column "source" lists the citations where the parameterizations were first defined. 



Force 


to 


ti 


t2 


t3 


Xo 


Xi 


X2 


x-s 


SkM* 


-2645.0 


410.0 


-135.0 


15595.0 


0.090 


0.0 


0.0 


0.0 


z,. 


-1983.76 


362.25 


-104.27 


11861.4 


1.1717 


0.0 


0.0 


1.7620 


SkT6 


-1794.2 


294.0 


-294.0 


12817.0 


0.392 


-0.5 


-0.5 


0.5 


bijy4 


o /I oo mo 
-2450. 9io 


480. OiO 


c /I f; one 


L6l 1 ( .0 


n oo /I n 
0.834U 


-0.3438 


1 A 

-i.U 


1 O K /I O 


Skll 


-1913.619 


439.809 


2697.594 


10592.267 


-0.954536 


-5.782388 


-1.287379 


-1.561421 


SkI3 


-1762.88 


561.608 


-227.090 


8106.2 


0.3083 


-1.1722 


-1.0907 


1.2926 


SkI4 


-1855.827 


473.829 


1006.855 


9703.607 


0.405082 


-2.889148 


-1.325150 


1.145203 


SkP 


-2931.70 


320.618 


-337.409 


18708.96 


0.29215 


0.65318 


-0.53732 


0.18103 


SkO 


-2103.653 


303.352 


791.674 


13553.252 


-0.210701 


-2.810752 


-1.461595 


-0.429881 


SkO' 


-2099.419 


301.531 


154.781 


13526.464 


-0.029503 


-1.325732 


-2.323439 


-0.147404 



Force 


64 




a 


h^/2m 




dcM 


Vp 




sc 


)U] 


-ce 


SkM* 


65.0 


65.0 


1/6 


20.7525 








279.082 


258.962 


I 


34 






61.845 


61.845 


1/4 


20.7525 


1 


1 


231.823 


222.369 


I 


36 




SkT6 


53.5 


53.5 


1/3 


20.750 


1 





202.526 


204.977 


I 


35 




SLy4 


61.5 


61.5 


1/6 


20.73553 








295.369 


286.669 


I 


32 




Skll 


62.130 


62.130 


1/4 


20.7525 





1 


285.209 


291.384 


I 


30 




SkI3 


94.254 


0.0 


1/4 


20.7525 





1 


335.432 


331.600 


I 


30 




SkI4 


183.097 


-180.351 


1/4 


20.7525 





1 


286.029 


310.832 


I 


30 




SkP 


50.0 


50.0 


1/6 


20.73 


1 





252.619 


236.237 


I 


37 




SkO 


176.578 


-198.7490 


1/4 


20.73553 





1 


253.771 


269.942 


I 


38 




SkO' 


143.895 


-82.8888 


1/4 


20.73553 


1 


1 


256.095 


258.122 


I 


38 





TABLE II. SMMC values of (Q^), {Ql), (Q^), and (QpQ„) 
(in 6"*) for 28,30,32,34 jyjg Typical error bar is ±26*. 







sdpf 






Nucleus 


iQ') 


{Ql) 


{Ql) 


(QpQn) 


^*^Mg 


29.6 


12.7 


5.0 


5.9 


3°Mg 


36.7 


14.5 


9.6 


6.3 


32 Mg 


42.6 


15.0 


14.5 


6.6 


3*Mg 


38.1 


12.2 


12.8 


6.5 






sdpf 






2«Mg 


37.2 


13.8 


11.7 


5.9 




30.0 


13.3 


8.7 


4.0 


32 Mg 


42.8 


13.5 


18.2 


5.5 


'"Mg 


68.1 


13.4 


35.7 


9.5 
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FIG. 1. The surface energy coefficient as defined in Eq. (B.13) for the Skyrme parameterizations used in this work. 
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FIG. 2. Potential energy surfaces for ' ' Mg as functions of quadrupole deformation /3 for the set of Skyrme parameteri- 
zations. 
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FIG. 3. Same as in Fig. |, except for 26,28,30,32 j^^^ 
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FIG. 5. Top: Spherical neutron shell structure in '^^Mg calculated in several Skyrme-HF models. The single-particle levels 
are normalized to the energy of the 0(3/2 shell. Bottom: The prolate-spherical energy difference calculated with (dashed line) 
and without (solid line) pairing. 
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FIG. 6. Top: Single-neutron dz/i and /7/2 energies as functions of Z predicted in several Skyrme-HF models. The posi- 
tive-energy levels are marked by a dashed line. Bottom: the size of the corresponding single-particle Ar=20 gap, Ae2o. 
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FIG. 7. Same as in Fig. |, except for 38>4o,42,44g^ 
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FIG. 8. Spherical neutron shell structure in S calculated in several Skyrme-HF models. The single-particle levels are 
normalized to the energy of the /7/2 shell (a); the position of the deformed prolate minimum with respect to the spherical HF 
state (b); the position of the deformed prolate minimum with respect to the spherical HF state (c). 
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FIG. 9. Same as in Fig. |, except for 80,82,842^. 
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FIG. 10. Spherical neutron shell structure in Zr calculated in several Skyrme-HF models. The single-particle levels are 
normalized to the energy of the pi/2 shell (a); spherical proton shell structure (b); the position of the deformed prolate and 
oblate minima with respect to the spherical HF state (no pairing) (c). 
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FIG. 11. Same as in Fig. |, except for 92.94,96.98,1002^. 
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FIG. 12. Spherical neutron shell structure in Zr calculated in several Skyrme-HF models. The single-particle levels are 

normalized to the energy of the pi/2 shell (a); spherical proton shell structure, the single-particle levels are normalized to the 
energy of the ^5/2 shell (b); the position of the deformed prolate and oblate minima with respect to the spherical HF state (no 
pairing) (c) 
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FIG. 13. The effect of ZPE corrections on the PES of ^^Mg calculated with the Skl4 SKyrme-HF model. Top two: proton 
and neutron pairing energies. Middle: uncorrected (dashed line) and ZPE-corrected (solid line) PES. Bottom: the various 

contributions to the ZPE (rotational, dashed lino; /9-vibrationaI, long dashed line; ccntcr-of-rnass, short dashed line) and the 
sum of rotational and vibrational corrections, solid line. The vertical lines mark borders of the regions where static pairing (in 
neutrons or protons) vanishes. 



29 




30 




-0.2 D 0.2 0.4 



FIG. 15. Same as in Fig. |l|, except for *°Zr. 
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FIG. 17. Results of SMMC calculations with the sdpf (left) and sdpf (right) interactions for the neutron-rich Ne and Mg 
isotopes: Top: (QpQn); Bottom: (A+A) (J=0, T=l). 
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FIG. 18. Single-neutron shell-model occupations in neutron-rich Ne and Mg isotopes calculated in the SMMC with sdp f 
(left) and sdpf (right) effective interactions. 
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FIG. 19. J=0, T=l matrix elements {jaja\V\j/3j/3} of the sdpf (open circles) and sdpf (dots) interactions. The diagonal 

matrix elements (q=/3) are indicated. The matrix elements are presented according to the following convention. Assuming 
that the single-particle orbitals are labeled as: ^3/2, •Si/2, ff/2, /5/2, P3/2, P1/2) = (1, 2, 3, 4 ,5, 6, 7), the first six matrix 

elements are those within the sd shell [(1,1), (1,2), (1,3), (2,2), (2,3), (3,3)], the next ten are the matrix elements within the 
pf shell [(4,4),(4,5)...,(7,7)], and the remaining twelve are the cross-shell matrix elements [(1,4), (1,5),...,(3,7)]. 
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FIG. 20. Single-particle ldg^2 and I/7/2 neutron levels ( B13 ) predicted in HFB-SM with the sdpf interaction for the Ne and 
Si isotopes with 14<Af<22. Note that the size of the N—20 gap changes very little with A*'. 
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FIG. 21. Single-particle energy gap at A'^=20, Ae2o, calculated in the HFB-SM with sdpf and sdpf' interactions for the 
A''=20 isotones of O, Ne, Mg, Si, and S. Note that in sdpf the size of the N—20 gap depends dramatically on Z. 
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FIG. 22. Comparison between the particle-hole matrix elements V{a,P), Eq. (B9), of sdpf and sdpf' . The same convention 
is used as in Fig. 

iTc=-Tp). 



19[ Top: neutron-neutron (proton-proton) matrix elements (ra—rp). Bottom: proton-neutron matrix elements 



38 




FIG. 23. Correlation energy (Mj) for the neutron-rich Ne and Mg isotopes. Top: sdpf results. Bottom sdpf' results. 
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FIG. 24. Size of the N—20 gap, Ae2o, obtained from canonical single-particle energies (dashed lines) and the shell-gap 
parameter obtained from two-neutron separation energies (solid lines) calculated in the HFB approach with Skyrme 

interactions SkP (dots) and SLy4 (circles). 
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